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Abstract 

We conduct a general discussion of supersymmetric models with three families in the Higgs 
sector. We analyse the scalar potential, and investigate the minima conditions, deriving the mass 
matrices for the scalar, pseudoscalar and charged states. Depending on the Yukawa couplings 
and the Higgs spectrum, the model might allow the occurrence of potentially dangerous flavour 
changing neutral currents at the tree-level. We compute model-independent contributions for 
several observables, and as an example we apply this general analysis to a specific model of 
quark-Higgs interactions, discussing how compatibility with current experimental data constrains 
the Higgs sector. 


1 Introduction 

Even though the standard model (SM) has proved to offer a very successful description of strong and 
electroweak interactions, it fails in providing an explanation to issues such as the gauge group, the 
number of families, the dynamics of flavour and the mechanism of mass generation, among others. A 
particularly puzzling topic is that of the number of quark and lepton generations. From a phenomeno¬ 
logical and experimental point of view, we have strong reasons to accept that there are indeed only 
three copies of up- and down-type quarks, as well of charged leptons and active neutrinos. However, on 
the theoretical side, there is little hint to the origin of this three-fold replication. The cancellation of 
anomalies requires identical number of quark and lepton families, but adds no information to the total 
number. Extensions of the SM, such as supersymmetry (SUSY), or grand unified theories (GUT), 
remedy several deficiencies of the SM, but they also fail in explaining the number of fermion families. 

Given the existence of three families of quarks and leptons, and since neither theory nor experiment 
impose any constraints on the number of Higgs families, one could wonder whether the same family 
replication takes place in the Higgs sector. In addition to being very aesthetic, this scenario is not 
unexpected from a theoretical point of view. In some SUSY models from strings, as is the case of the 
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compactification of the 10 dimensional Eg x E' s heterotic superstring to 4 dimensions, the low-energy 
theory thus obtained is associated with a subgroup of the gauge group Eg. Under Eg, fermions and 
Higgses are assigned to the same representation, and the requirement of three families of fermions 
imposes replication of the Higgs content [1]. In fact, many string constructions that contain three 
fermion families also include Higgs family replication, as is the case of Zg orbifolcl compactifications of 
the heterotic superstring, which also exhibit a low-energy spectrum composed of three supersymmetric 
Higgs families [2-7], and that of some D-brane models [8]. Other models that include a non-minimal 
Higgs content have also been extensively addressed in the literature [9-25]. The consequences of 
extending the Higgs sector are abundant, and have implications which range from the theoretical to 
the experimental level. For instance, and if the extra Higgses are light, the addition of these states in a 
minimal SUSY scenario will spoil the unification of gauge couplings around 10 16 GeV. Nevertheless, in 
models from the heterotic string, the high energy scale is different (« 5gut x 10 17 GeV), and the new 
states can even be helpful regarding unification [5]. Since ours is a low-energy oriented analysis, we will 
not consider the (more speculative) high-energy implications of the extended Higgs sector. The most 
challenging implication of an extended Higgs sector is perhaps the potential occurrence of tree-level 
flavour changing neutral currents (FCNC), mediated by the exchange of neutral Higgs states 1 . In the 
SM and the minimal supersymmetric standard model (MSSM), these effects are absent at the tree- 
level, since the coupling of the quark-quark-Higgs mass eigenstates is flavour conserving. This arises 
from having the Yukawa couplings proportional to the quark mass matrices, so that diagonalising 
the mass matrices also diagonalises the Yukawas. Since experimental data is in good agreement with 
the SM predictions, the potentially large contributions arising from the tree-level interactions must be 
suppressed in order to have a model which is experimentally viable. In general, the most stringent limit 
on the flavour-changing processes emerges from the small value of the K^ — Ks mass difference [10,11]. 

Avoiding the FCNCs induced by the tree-level exchange of neutral Higgs bosons can be achieved 
by several distinct approaches, each involving a different sector of the model. 

(i) Discrete symmetries. Imposing a discrete symmetry on the model ensures that only one generation 
of Higgs couples to quark and leptons, and completely eliminates tree level FCNCs from the predictions 
of the model [29,30]. Naturally, one attempts to motivate such an assumption via topological and/or 
geometrical arguments. For instance, in [16,17], adding four extra Higgs doublets to the MSSM 
content, and assuming that the FCNCs thus induced are suppressed via some symmetry, then the 
extra generations (labelled “pseudo-Higgs bosons”) do not acquire vacuum expectation values (VEVs), 
do not mix with the two MSSM-like doublets, nor couple to fermions. Still, the lightest of the new 
states is stable so that it becomes a candidate for dark matter. 

(ii) Suppression of Yukawa couplings. In this case, FCNCs are not eliminated, but the Yukawa 
interaction responsible for the tree-level flavour violation is made small enough to render the new 
contributions negligible. Since the most stringent bounds are usually associated with tree level con¬ 
tributions to the neutral kaon mass difference (A mu), the Yukawa couplings of down and strange 
quarks are forced to be very small. For example, in [25], a U( 2) flavour symmetry was considered as 
the candidate symmetry to suppress FCNCs. Although this second possibility is more appealing in 
the sense that it does not require excluding part of the Higgs content from Higgs-matter interactions, 

1 For other consequences such as flavour-violating Higgs and top decays, and the associated experimental signatures 
at the next generation of colliders, see, for example, [26—28] and references therein. 
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its major shortcoming lies in the fact that in general one lacks a full theory of flavour, which would 
predict the Yukawa couplings 2 , as would be the case of string theory. 

(iii) Decoupling of extra Higgses. If one does not wish (or is not allowed) to impose a symmetry, or 
if the Yukawa couplings are not free parameters of the model, a third possibility lies in making the 
new Higgs states heavy enough, so that the contributions they mediate are suppressed [9-11,13,19]. 
The new states are thus decoupled, and the effective low-energy theory is very similar to the usual 
MSSM/SM. This possibility is sometimes the only “degree of freedom” remaining, especially in highly 
predictive models, where the Yukawas arise from some high-energy formulation, as is the case of string 
models. Nevertheless, it is important to stress that the decoupling scenario is in general achieved by 
enforcing very large values for some of the Higgs soft-breaking masses. This may lead to a fine-tuning 
scenario in association with electroweak (EW) symmetry breaking. When the decoupling approach is 
taken together with the suppression of Yukawa couplings (ii), it is possible to obtain Higgs spectra 
which manage to comply with experiment without excessively heavy Higgses [10,11,13,19-21,25]. 

From the above discussion, it is clear that a correct evaluation of the FCNC problem associated to 
multi-Higgs doublet models is indeed instrumental. In what follows, we propose to analyse the most 
general scenario of the MSSM extended to include three families of Higgs doublets. We do not impose 
any symmetry on the model, allowing for the most general formulation of both the superpotential 
and the SUSY soft breaking Lagrangian. We provide a generic overview of the extended Higgs sector, 
and discuss the minimisation of the scalar potential. In this work, we also compute the most general 
expression for the contribution of tree-level neutral Higgs exchange to neutral meson mass difference. 
Contrary to previous analyses [10,11,13,19], we include the contributions from all physical (rather 
than interaction) Higgs states, scalar and pseudoscalar. Moreover, we take into account the mixing 
in the Higgs sector. Finally, as an example of how to apply our general formulation, we consider an 
ansatz for the Yukawa couplings along the lines of previous analyses (namely the “simple Fritzsch 
scheme” [31-33]), and evaluate the specific contributions to the neutral mesons mass differences. Our 
findings turn out to be more severe than those of previous works. 

The work is organised as follows. In Section 2, we analyse the extended Higgs sector, paying special 
attention to the minimisation of the Higgs potential and addressing potential fine-tuning issues. We 
compute the tree level mass matrices in Section 3 and discuss the associated spectra. Section 4 is 
devoted to the analysis of Higgs-matter interactions, and in Section 5 we derive a model independent 
computation of the tree-level contributions to the neutral meson mass differences. Assuming an 
illustrative example for the Yukawa couplings, we present in Section 6 a numerical study of the FCNC 
contributions, investigating how the Higgs spectrum must be constrained in order to have compatibility 
with experimental data. Finally we summarise the most relevant points in Section 7. 


2 Extended Higgs sector 

We begin our analysis by addressing the Higgs sector of a SUSY model where three generations of 
SU( 2) doublet superfields are comprised. In each generation, one finds hypercharge —1/2 and +1/2 

Nevertheless, it is worth noticing that in the framework of the MSSM, there are several models where it is found 
that via flavour symmetries one can successfully explain the observed pattern of masses and mixings. 


3 



fields, coupling to down- and up-quarks, respectively: 


H 


1(3,5) 


^1(3,5) 

^1(3,5) 


H 


2(4,6) 


""2(4,6) 

h° 

"2(4,6) 


2.1 Tree-level Higgs potential 


(1) 


The most general superpotential of a model with three families of Higgs doublets can be written as 
follows: 


W = Q (Y d H t + Y d H 3 + Y d H 5 )D c + L {Y^H X + Y 3 H 3 + Y£H 5 )E C 
+ Q {Y 2 u H 2 + Y^Hi + Y£H 6 )U c + i-iv2H[H'2 + mHiH 4 + Mie^i H 6 
+ M 32 H 3 H 2 + / 134 H 3 H 4 + h 3 qH 3 H§ + n§ 2 H 3 H 2 + (£ 54 H 3 H 4 + /X 56 H 3 Hq . (2) 

Above, Q and L denote the quark and lepton SU(2)l doublet superfields, U c and D c are quark 
singlets, E c the lepton singlet, and Y® are the Yukawa matrices associated with each Higgs superfield. 
Regarding the /x-terms, these are now extended in order to include all possible bilinear Higgs terms. 
Calling upon some specific model that would be responsible for generating the latter terms, it would 
be possible to reduce the number of free parameters. It is also important to recall that one could also 
impose a (discrete) symmetry acting on the superpotential, whose effect would be the suppression of 
some of the couplings in W. However, in our discussion, we consider the most generic form for W, 
thus taking the several bilinear parameters in W as effective /u-terms. 

From the above one can derive the F- and D- terms, writing the latter in doublet component for 
simplicity: 


Vf = Y, VuNiH}H j+ Y tiHikHtm, 

i,j— 1,3,5 £=1,3,5 

1= 2,4,6 k,l—2,4,6 


J2 3 1 

" 6 

2 

„/2 

" 6 

wi 

^i°° 

H 

£ 

Y h I t ° h " 
.£=1 

+ V 

Ec-imi 2 

_£=1 


where Hi are the scalar doublets belonging to the superfields of Eq. (1), g, g' denote the gauge coupling 
constants, and r“ are the SU(2)l generators. The soft SUSY-breaking terms trivially generalise the 
minimal two-doublet case, including the usual soft-breaking masses and Bfi -like terms in the Higgs 
potential. 

y soft = Y H \ Hj + Y (mDkiHtHi- Y [{B fi ) il H i H l +K.c.]. (4) 

£,.7=1,3,5 k,l— 2,4,6 £=1,3,5 

l—2,4,6 

2.2 Minimisation of the potential: the superpotential basis and the Higgs 
basis 

After electroweak symmetry breaking, the neutral components of the six Higgs doublets develop the 
following VEVs, 

(*1(3,5)) = W 1(3,5J j (*2(4,6)) = ^2(4,6) i (5) 
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and as usual one can write 


ti 


uu 


i , 

7i (OT 


■ iVi) ■ 


( 6 ) 


In the above we assume that all the VEVs are real 3 . The next step is to minimise the scalar potential 
with respect to the VEVs. Although the problem appears to be a simple generalisation of the two- 
Higgs doublet model (MSSM), it will prove to be much more involved. It happens that in the presence 
of six non-vanishing VEVs, finding a solution to the minima equations is a rather cumbersome task. 
As an illustration, we present the equations obtained from minimising with respect to the down-type 
Higgs: 


E 


( m rf)A+ E f^jl ftil 

Z=2,4,6 


1= 2,4,6 


wi + wj i Mf 


= 0, 


where j = 1,3, 5 and we have defined 


(7) 


Vd, = 



Vu 



( 8 ) 


In addition to Eq. (7), which only ensure that we are in the presence of an extremum, we have 
to further impose the conditions for a minimum with respect to the several variables. Therefore, 
and for the down-type Higgs fields, having positive second derivatives is equivalent to the following 
inequalities: 


(' m d)jj + E Mj£ + 9 


v 2 , — vJ.+2 w 2 - 


> 0, 


(9) 


£=2,4,6 


provided that at the minimum the determinant of the six-dimensional matrix is positive, or equiva¬ 
lently, det(<9 2 V / Oh® dh *•) > 0. 

The equations for the remaining cases (minimising with respect to the up-type Higgs fields) can 
be obtained from Eqs. (7,9) by interchanging {1,3,5} <-> {2,4,6} and changing the sign in the terms 
proportional to Mf as Mf (v d — v 2 ) —» —M\[v 2 d — v 2 ) (cf. Eq. (7)) and Mf (v d — v 2 + 2 w 2 ) —» 
—Mf (v d — v 2 — 2w 2 ) (cf. Eq. (9)). The only choice of parameters that would render the solution of 
the minima equations straightforward is to call upon the diagonal soft breaking masses (( md)u and 
( m u )kk ) to be fixed by the minima equations. For (rrid)\\ the minimum equation would read 


K)u = E — - ( TO d)i3 — - ( m d )is — 

. “« w i w i 

3- 2,4,6 


W 5 


- E 


(E 


*=1,3,5 Z=2,4,6 


... .. W M l * / rl v d- v u 

Ml lVil - ) — 9 M Z -2 

Wx 2 V 2 d + V 2 


( 10 ) 


Naturally, the equations can be solved for a distinct set of parameters, but then finding an analytical 
solution is in general not possible. Analytical relations are very useful when writing down the mass 
matrices for the several Higgs sectors, since one can directly replace the parameters. In fact, the 

3 We do not discuss the possibility of spontaneous CP violation in association with this class of multi-Higgs doublet 
models. 
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major handicap when relying on the diagonal soft masses as the minimum parameters is that then 
one loses control of the leading contributions to the Higgs boson mass eigenvalues. This will be made 
clear in a forthcoming section. Another shortcoming is that, even though we can define a generalised 
expression for tan (3 as tan /3 = \Jw\ + w 2 + Wq / \Jw\ + iff + , writing the minimisation equations 

as a function of sin/3, cos/3 (in an MSSM-like fashion), is impossible. 

So far we have been working in the basis that naturally emerges from both the superpotential and 
soft-breaking Lagrangian formulation. To simplify the discussion, we will henceforth label this natural 
basis as the superpotential basis. However, this need not be the unique approach when addressing 
the minimisation of the potential. 


As it has been shown, one can work in a basis where only two of the new six neutral fields have 
non-vanishing VEVs, the so-called Higgs basis [9,15]. For the down-type Higgses, the new fields, fa, 
are related to the original ones by means of the following unitary transformation, P d '. 




/ 9 

V w l 


Wo 




(wi h\ + u >3 hi + w 5 h°) 


( 11 ) 


b° 3 = 


/ 9 

V w i 


Wo 


Wr 


Wo + Wr hi — 


W 1 


y/ w 3 + w l 


(w 3 hi + w 5 h° 5 ) 


( 12 ) 


</>5 = 7 o 1 ( w 5 h° 3 - W 3 h° 5 ) . 

V w 3 + w 5 

An analogous transformation, P u , can be derived for the up-cjuark-coupling Higgses 
adequate replacements ({1,3,5} —> {2,4,6}). For simplicity, let us introduce the 
parametrisation for the transformations Pd and P u : 


(13) 

{hi 4 6 ), with the 
following global 


— P h 


(14) 


where P t j is a 6 x 6 matrix whose entries are defined from Pd and P u as 


P = 

1 13 


{Pd)i+±i±L i,j = 1,3,5 

{ p u)^L = 2,4,6 

0 elsewhere. 


From the above, it is clear that in the new basis only two of the fields do have a VEV 

(</>?) = v d , ( 0 °) = v u , 

with Vd, u defined in Eq. (8), which in turn must satisfy 

vl + v 2 d = 2 M 2 z /{g 2 + g' 2 ) « (174 GeV) 2 , 


(15) 


(16) 


(17) 


where Mz is the mass of the Z° boson. In this basis one has the additional advantage that, similar 
to what occurs in the MSSM case, one can clearly define tan /3: 


We can now write 


tan (3 = — 
Vd 


1 


(18) 


^1,2 v d,u + ~^= {R\,2 + *^ 1 , 2 ) 


i = 3-6. 


(19) 
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Thus, one can now understand the three Higgs family model as an MSSM-like model, extended by 
4 additional doublets, which do not directly interfere with the mechanism of electroweak symmetry 
breaking. In the Higgs basis, the potential preserves its original structure with respect to the depen¬ 
dence on the field combination, but the parameters associated with the F- and soft breaking terms 
are now redefined as 


where 


V F = ^2 (Ab <I>1 fa + Mcd< 


a,b= 1,3,5 
1=2,4,6 


i= 1,3,5 
c,d= 2,4,6 


Kioft — ^ ^ ( V^d)cib fi\ fib + ^ ^ ( V^u)cd fi\ fa- Y. ac fia fic + H.c.] 


a,b= 1,3,5 


c,d= 2,4,6 


a=l,3,5 
c=2,4,6 


Mail = ( Pd)ai l^il Hjl i^d)jb , 
(jTld)ab = ( Pd)ai ( m d)ij (-^d)i6 > 
(3?/x) ac = (pj)ia ( Pu)lc , 


Hcd (Hu)ci Mii Mifc {Pu)kd , 

(^u)cd = ( Pu)ck ( m u)kl {Pu)ld > 

i,j,a,b= 1,3,5; k,l,c,d= 2,4,6. 


In the new basis, the problem of minimising the scalar potential is quite simplified, and one is able 
to derive the six minimisation conditions in a compact form. The conditions for the down-type sector 
read 

2 M 7 

mu = bi 2 tan /3- — cos 2/3 , 

m 13 = ^32 tan /3 , 

m i 5 = ^52 tan/3, (22) 

while those related to the up-type Higgses are given by 

M 2 

772-22 = &12 COt /3 + —COS 2/3 , 
m 24 = bu cot /3 , 

m 26 = &16 c °t /3 • (23) 


In the above, and for simplicity, we have introduced the short-hand notation 


m ij 


E/= 2,4,6 My- + (m 2 d)ij i, j — 1,3,5 
Efe=i,3,5 My- + ("ij)y i,j = 2,4,6, 


bij = (HyLt)j 


where we stress that both parameters /Xy, and /Xy, as well as fry have dimensions (mass 2 ). Throughout 
the paper we will be often using the above short-hand notation, even though this implies that the 
rotated /x-terms and Higgs soft breaking masses loose their individual character, becoming merged 
into the quantities m 2 ; -. As a final comment, let us notice that the minima equations (22,23) are, in 
structure, very similar to those of the MSSM, namely the first in each set of three. It is worth referring 
that the EW scale is only explicitly present in two equations (those derived with respect to </>i 2 , the 
VEV-acquiring fields). Notice however, that this is a mere consequence of working in a different basis. 
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In addition to the minima conditions Eqs. (22,23), one can also derive other conditions, which are 
useful in constraining the allowed parameter space. As an example, and from direct comparison of 
the minima equations for mf 1 and m| 2 > we find the following inequality: 

m ll™22 < b 12, (25) 

strongly resembling the analogous MSSM condition. 

2.3 Numerical minimisation of the potential: illustrative examples 

In what follows, our goal is to conduct a short analysis of the issues discussed in the previous sub¬ 
sections. The phenomenological analysis of the Higgs sector (masses and mixings) will be carried in 
detail in Section 3, working in the Higgs basis. Nevertheless, we find it interesting to evaluate how the 
minimisation of the scalar potential constrains the soft breaking terms. To do so, we consider several 
numerical examples for the minimisation of the scalar potential for distinct VEV regimes. In addition, 
we also discuss whether or not a specific choice of soft breaking terms may translate into a fine-tuning 
problem. Although the minima conditions are cast in a much more compact and appealing way in 
the Higgs basis, since only two of them explicitly involve the EW scale, a possible fine tuning problem 
may be unapparent in this case. Therefore, for this specific discussion, we will stick to the original 
superpotential basis. In this case, and since we make no assumptions regarding the mechanism of 
SUSY breaking associated with this generic multi-Higgs doublet model, we take the input parameters 

( m d)ij > ( m u)ij r L l ij > , (26) 

to be free at the electroweak scale. Regarding the Higgs VEVs (to,), we impose no constraint other 
than satisfying the EW breaking conditions for a given value of tan (3. In order to further simplify 
the approach, we will fix, in each separate analysis, the value of (ijj = /i, and assume a common 
overall scale for the off-diagonal masses ( rrid)ij and ( m u )ij , and for (. More specifically, we take 
universal Bfi terms, (R/i),, = R/i = M§, where Ms is a “typical soft-SUSY breaking scale”. For the 
off-diagonal masses ( m,d)ij and ( m u )ij , the common scale is taken to be ( m.d)ij = ( rn u )ij = 0.1 Mg 4 . 

As discussed in the previous subsection, the parameters which are (trivially) fixed by the six 
minima equations are the diagonal Higgs soft-breaking masses, (mfya and (mqju- Clearly, from 
Eq. (10), the minima values for (m^) a, and in particular those of the down-type Higgs sector, will 
be very dependent of the values of the VEVs considered, most specifically on whether one chooses a 
degenerate or hierarchical regime. Therefore, each case has to be separately investigated, and in each 
situation, one must analyse the impact of the other free parameters, namely that of the overall scale 
chosen for the soft-breaking parameters. 

In addition, one has to investigate whether or not the parameters taken lead to a potential fine 
tuning (FT) problem for the parameters. In doing so, our approach will be the following. For a 
specific scenario defined by the parameters in Eq. (26), together with Wi and tan/3, we compute the 
values of the diagonal soft breaking masses, as imposed from complying with the minima equations. 

4 Notice that for i y j, values of ( fn q )ij /// Mg are not compatible with viable solutions. This is evident from the 
minimisation conditions of Eq. (10), using > 0- Furthermore, the limit ( m q )ij —> 0 occurs in well motivated 

models [34]. 



We then impose a tiny perturbation on the solutions, 


K)ii -> {m 2 q y u = (1 + A 9 ) ( m 2 q )u , 


(27) 


where A is taken to be around a few percent. Finally, we compute the new value of Mz derived from 
perturbing the correct (true minima) solution for ( [nig)a , M' z . To evaluate the amount of FT, we 
follow [35,36], and introduce the parameter A p, defined as 



5 Ml 


(28) 


identifying pi with ( m 2 )u . A rough measure of the fine-tuning can be derived from (Apj) -1 , since the 
latter can be identified with the probability of a cancellation between terms of a given size to obtain 
a final result Api times smaller. 

In what follows we take, as illustrative examples, two distinct cases for the Higgs VEVs: degener¬ 
ate and hierarchical wy. We stress here that at this point we are only investigating the constraints on 
the soft breaking terms arising from the minima conditions. A study of the spectrum is far simplified 
when moving to Higgs basis, as will be done in the following Section. 


• Degenerate VEVs: wi = W 3 = w$, w 2 = W 4 = w& 

In this case, the minima equations become much simpler, since in each sector the VEVs factor out. 
Therefore, having fixed the free parameters as above described, the six original equations essentially 
reduce to two independent ones: one for the down-type Higgs (e.g. ( fnd)\i ) and another for the 

up-type sector (e.g. (m u ) 22 )- 

In Fig. 1, we plot {jn ( i)\i and (m u ) 2 2l as computed from the tree-level minima equations, as a 
function of the common soft breaking term scale, Ms, for p. tJ = p = 200 GeV. In each case, we 
consider two distinct values for tan/3: 5 and 10. From Fig. 1, several properties of this model become 
apparent. First of all, it becomes clear that in this basis (the superpotential basis) finding minima can 
be a challenging task. Notice that although one can find solutions for the down-type sector associated 
with low Ms, in this range we are in the presence of maxima (rather than minima) for the up-type soft 
masses. The situation slightly worsens with increasing tan /3. Taking other (higher) values of /i would 
lead to the displacement of the vertical line further to the right of the parameter space parametrised 
by Ms- This behaviour can be confirmed by inspection of Eqs. (7,9,10). 

Let us now evaluate how fine-tuned these solutions are. In Fig. 2 we investigate A p for the up 
and down sectors, by studying in each case the effect of perturbing the minima solutions (m ^) 2 and 
(m 22 ) 2 by \ d,u = 1% and A d,u = 5%. We again consider p = 200 GeV, and take tan/3 = 10. 

From Eqs. (27,28) together with the minimisation conditions in Eq. (10), it is apparent that 

( 2 V 

A pi oc Moreover, it can be clearly seen that for the down-type sector, the dominant term (Bp) 

is further enhanced by a factor of tan/3, so that a tiny fluctuation in (m ^) 2 is not easily cancelled. 
In the up-type masses, the situation is reversed. The up-type version of Eq. (10) would exhibit a 
suppression of Bp by cot /3, so that here we find a more relaxed scenario. 

• Hierarchical VEVs: w 5 = 10 w 3 = 100 W\, wq = IO 1 U 4 = 100 w 2 
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Figure 1: Minima solutions mf ± and as a function of the common SUSY breaking scale Ms for 
tan p =5, 10 and degenerate VEVs. To the left of the vertical dotted line, the solutions for m'f are 
not true minima of the potential. 



Figure 2: A p as a function of the SUSY scale Ms for perturbations of 1% and 5% in ( m f-J 2 and 
( m 22 ) 2 - A 4 = 200 GeV, and tan p = 10 (the two lines for Apn appear over-set). To the left of the 
vertical dotted line, the solutions for are not true minima of the potential. Horizontal lines denote 
areas above which one roughly expects a FT stronger than 1% and 10%. 


In this case, and when compared to the degenerate one, each soft mass exhibits a very distinct 
behaviour at the minima of the potential. We will thus analyse the solutions which in each case are 
associated with the lightest and heaviest VEV, {mo) n, (m u )| 2 and (md)^, (m u ) qq- As before, in 
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10 5 
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Figure 3: Minima solutions for non-degenerate VEVs as a function of the SUSY scale Ms, with 
/r = 200 GeV, and tan/3 =10. To the left of the vertical dotted line, the solutions for mf, are not true 
minima of the potential. 


Fig. 3 we display the above masses as a function of Ms for /r = 200 GeV, and tan/3 =10, similar to 
what was presented for degenerate VEVs in Fig. 1. As seen from Fig. 3, the range for the solutions is 
not strongly affected. Notice however that the minima values of (m u )f 2 and become degenerate 

for large values of Ms- Moreover, the soft mass associated with the largest VEV (mg 6 ) presents, as 
expected, much smaller values than the others. This is a consequence of the suppression of Wj/we on 
B[i. Regarding the fine tuning, the situation is more involved. In Fig. 4 we plot A p for the up and 
down sectors, considering only the effect of a X d ’ u = A = 1% perturbation in (m 22 ) 2 , (m| 4 ) 2 , (mf 5 ) 2 
and (mgg) 2 . We take n = 200 GeV and tan/3 = 10. 

As seen from Fig. 4, the alterations emerge in both up- and down-type Higgs sectors. Regarding 
the down-type scalar (A p^) the fine tuning is typically larger now, and the effect is essentially due 
to the fact that it is now much more difficult to find cancellations among the several terms entering 
in the minimisation equations. In the up-type sector, we find a very similar situation to that of Fig. 2 
for the second and third generations (associated with the larger VEVs). Conversely, and since it is 
now associated to a tiny VEV, the first generation up-type diagonal soft mass is more unstable under 
perturbations, and its behaviour merges with that of for larger values of Mg. 

It has been shown in the MSSM that with a soft SUSY scale of a few hundred GeV the associated 
FT is around the level of 10% [36]. When compared to the MSSM, this extended model offers a more 
problematic FT scenario, as is manifest in Figs. 2 and 4. This can also be seen from the comparison 
of Eq. (10) with its MSSM counterpart. In addition to the single (Bfi) 12 and /i 42 terms, one now 
encounters additional (B[i)ij and Hufiu, as well as new soft breaking masses on the right hand-side 
of the above equation, which must cancel out to achieve the correct value of Mz- 

Finally, and even though in our analysis the VEVs were not fit by the minima equations, but taken 
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Figure 4: A p as a function of the SUSY scale Ms for perturbations of 1% in the diagonal soft masses, 
for non-degenerate VEVS. We take p =200 GeV and tan/3 =10. To the left of the vertical dotted 
line, the solutions for rn/1 are not true minima of the potential. Horizontal lines denote areas above 
which one roughly expects a FT stronger than 1% and 10%. 


as input parameters instead, let us consider the effect of imposing a small perturbation on the VEVs, 
while keeping the other parameters as defined either by input values and/or minima conditions. As 
for the previous case of perturbing the minima values of m^, we now impose that 

wi —> w[ = (1 + p q ) Wi . (29) 

In Fig. 5, we plot the effect on Mz (parametrised by A p) arising from taking p d = p u = p = 1 % for 
the case of degenerate (left) and non-degenerate VEVs (right). Again we consider p, = 200 GeV and 
tan/3 = lO.igs/ 

It suffices to comment that the larger the VEVs, the more stable they are under perturbations, 
and this effect is particularly manifest for the non-degenerate VEV case. 


3 Tree-level Higgs mass matrices 

We are now in conditions to investigate the Higgs spectrum, which will contain eleven neutral states, 
and ten charged physical particles. Let us begin with the derivation of the tree level mass matrices 
for the charged scalars. 
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Figure 5: A p as a function of the common SUSY scale Ms for perturbations of 1% in Wi , for /r =200 
GeV and tan /3 = 10. On the left we present the degenerate VEV case, on the right the non-degenerate 
case. 


3.1 Charged Higgses 

In the basis defined by 5 (0(~, 02~*, 03~, 04*> 05 ~, 06 )i the scalar potential for the charged components 
reads, 

.2 6 


U±= ^ m?|0r| 2 + Y, m 2 i \4>t\ 2 + 9 —^Y\H i \ 2 \H i \ 2 + y - Y\ H l H i\ 


i=l,3,5 


i=2,4,6 


2 i /2 6 

9 + 9 


E \ H ^\ H i\ 2 - 


g 2 - 9 


i —1 

/2 6 


*>j=l 




*i3=l 

z-)-j=odd 


*J =1 

2+ji=even 


+ 0]"* 0^ + m| 5 03 * 0^ + m| 4 0^* 0| + 02 * 06~ + ™46 04 * 06 

+ H.C.) + (&12 0r 02~ + *>14 0r 04 + & 16 01" 06 + *>32 03 0l + *>34 0J 0 4 

+ *>36 03~ 06 + *>52 05 0l + *>54 0^ 04 + *>56 0 5 06 + H ' C ') ■ 

At the minimum of the potential, the mass matrix for the charged states is: 


Ml = 


^ *>12 tan 0-\-Myy sin 2 (3 
b 12 + ^ Myy s * n 2/3 

*> 32 tan /3 

*>14 

*>52 tan /3 

V b 16 


b 12 + ^ M vU sin2 ^ 
*>12 cot /3-\-Myy cos 2 (3 
b 32 

*>14 c °t /3 
*>52 

*>16 cot (3 


*>32 tan (3 
b 32 

m 33 ~T\V 

*>34 

m 2 
m 35 

b 36 


°14 

*>14 cot (3 
b 34 

m 44 + T ’W 
b 54 

m 2 
m 46 


*>52 tan (3 

b 52 

2 

m 35 

*>54 

2 rji 

m 55 ~ T W 
b 56 


*>16 ^ 
*>16 cot 0 
b 36 
m 46 
b 56 

m 66 + T W / 


(30) 


(31) 


where Miy is the mass of the W ± boson and TV = cos(2/3)(l — tan 2 6 w)- One can easily identify 

the massless combinations that will give rise to the charged Goldstone bosons (G 1 * 1 ), “eaten away” as 

5 The rotated basis of the charged fields is obtained via an identical transformation of Eq. (14). 
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the W ± acquire mass. The eigenstates will be henceforth denoted by hf, with i = 1 corresponding to 
the unphysical massless state. The upper-left 4x4 sub-matrix of the one above displayed is similar 
in structure to that derived in [15], in the framework of a four-Higgs doublet model. 

We stress here that the compact appearance of the matrix in Eq. (31) is a direct consequence of 
working in the Higgs basis. From direct inspection of the above matrix, one can impose that the 
diagonal blocks have semi-positive eigenvalues, thus deriving a zeroth order condition for avoiding 
tachyonic charged states. For example, for the {<j >3 — <f >4 ) sector, this would read: 

TO 33 m 44 £ b 34 , (32) 

while an identical argumentation within a sector (e.g.(0[j~ — <j>§) in the down-type Higgs) would in 
turn lead to 

m 33 to 55 > m 35 • (33) 

Albeit very useful, the above equations are only necessary (rather than sufficient) conditions to ensure 
the presence of true minima, so they should be interpreted as a means of orientation in the parameter 
space. Additionally, we notice that simple conditions as those above are only possible to derive in the 
Higgs basis. 


3.2 Neutral Higgses 


In our work, we are particularly interested in the neutral Higgs bosons. In order to find the neutral 
spectrum of the model, we decompose the complex fields as in Eq. (19), and derive the mass matrices 
associated with the real and imaginary components. CP conservation in the Higgs sector translates 
in the absence of terms mixing the Rt and I t components, so that scalar and pseudoscalar states do 
not mix. For the scalar mass matrix we obtain 





' H 2 dRidRi 









" 1 

min 




/ 

fc>12 tan (3-\-Mg cos^ (3 

-*>12 - 2 sin 2/3 

*>32 tan P 


-bl4 

*>52 tan P 

-*>16 

\ 


Mli 

-*>12 - sin 2(3 

&12 cot /3+M^ sin^ (3 

-*>32 

1 

cot (3 

-*>52 

*>16 co * P 



*>32 tan (3 

-*>32 

0 

TTT-gg-l - 2 — COS 2/3 


~ b 34 

m 35 

-*>36 



— *>14 

* > 14 co * ; P 

~ b 34 

™2 

m 44 

- ^2- cos 2(3 

-*>54 

m 46 



*>52 tan P 

~ b 52 

m 35 


-*>54 

9 

m 55 + ~T~ cos 

“ b 56 


V 

-*>16 

i>16 cot /3 

~ h 36 


m 46 

“*>56 

9 M % 

m 66 - T~ cos 2/3 

/ 


(34) 

It is straightforward to recognise the MSSM scalar Higgs mass matrix as the 2x2 upper left block of 
the previous 6 x 6 matrix. Likewise, the tree-level mass matrix for the pseudoscalar can be written as 


M 


2 

1 


1 d 1 2 V 


2 eiidij 


min 


14 



i >12 tan (3 
£>12 

£>32 tan P 
£>14 

£>52 tan P 
£>16 


°12 

£>12 cot P 

£>32 

£>14 co * P 
£>52 

£>16 cot P 


£>32 tan P 
£>32 

n 33^~ ! T~ cos 2 P 
£>34 


°14 

£>14 c °t P 
£>34 

2 M 7 

»'44 <2 cos 2 P 

£>54 

2 

m 46 


£>52 tan P 

£>52 

2 

m 35 

£>54 


5 + ^ cos 2/3 
£>56 


£>16 cot P 

£>36 
2 

m 46 
£>56 

M 7 / 

6 - 2 ^" cos 2/3 / 


( 35 ) 


In the 2x2 submatrix defined by the i,j = 1,2 entries, it is easy to identify the combination associated 
with the massless Goldstone boson. As in the MSSM, the latter degree of freedom is “eaten” by the 
Z° boson, as it acquires a mass. After rotating away the massless state, 

G° = —-j= (h sin j3 — I\ cos (3 ), (36) 

v 2 

by means of the unitary transformation 


( 7 r \ _ ( cos(3 -sin/3 A / h \ 

Y I' 2 J y sin/3 cos/3 JyhJ 

the five remaining eigenstates are in general massive, and the new mass matrix reads 


/ 

_2£>12 
sin 2/3 

b 32 sec /3 

b 14 esc (3 

652 sec (3 

bie esc (3 

\ 


632 sec (3 

777.2 

m 33“l- * T ~ COS ^ 

634 

m 35 

&36 



614 esc (3 

&34 

777, 2 

m 44-2^" COS ^@ 

&54 

m 46 



652 sec (3 

m 35 

&54 

m 2 

m 55“l— T ~ cos 2/3 

£>56 


V 

bie esc /3 

t>36 

m 46 

&56 

m 2 

m 66- f- COS 

2/3 J 


(37) 


(38) 


When discussing the issue of minimising the scalar potential, we pointed out that each of the presented 
basis (Higgs and superpotential basis) had its own advantages/drawbacks. Here we are facing a very 
strong point in favour of the Higgs basis, namely the possibility of rotating away the pseudoscalar 
massless state via a transformation that only involves the first two states. Should we have been working 
in the superpotential basis, the transformation of Eq. (37) would be far more complex. Even though 
one could still define tan /3 as tan /3 = \Jw\ + w\ + Wq / \Jw\ + ref + w 2 , the rotation parametrised 
by /3 would be extended to a 6 x 6 matrix, which would act upon the combination of Pi (/ 1 , 12 , ■ ■. Ie) T ■ 


It is worth recalling that the above mass matrices (scalar and pseudoscalar) are not those associated 
with the original (interaction) eigenstates. The relation between the superpotential basis and the Higgs 
basis is trivially obtained from Eq. (14), 


Ri — PijGj , I'i — Pjjj , i,j — 1 ... 6 . (39) 

In the Higgs basis, the mass matrices can be easily diagonalised by 

SrMrSI, = A 2 r = diag (mf ), 

SiM]s\ = A'j = diag(?nf 2 ), (40) 

where 7 are the diagonal scalar and pseudoscalar squared mass eigenvalues (notice that the i = 1 
term for the pseudoscalars corresponds to the unphysical massless would-be Goldstone boson). From 
the above it is straightforward to observe that the matrices that diagonalise the mass matrices in the 
original basis can be related to the latter as 


Sa,<p = Sr j P. 


(41) 
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3.3 Tree-level Higgs spectrum - a brief discussion 


Although a thorough analysis of the Higgs parameter space lies beyond the scope of this paper, it 
is important to comment on a few issues. The first regards radiative corrections, which play a key 
role in the MSSM and its extensions, especially in relation with the mass of the lightest scalar Higgs. 
Without radiative corrections, and as occurs in the MSSM, the mass of the lightest scalar is bounded 
to be rri h o < Mz [12]. In the framework of supersymmetric models with extended Higgs sectors, 
higher order contributions to the lightest Higgs mass have been already analysed [37], and the upper 
bound on m h o does not differ significantly from the one derived in the MSSM. Since in this work our 
aim is not so much to study in depth the Higgs spectra, but rather investigate to which extent FCNCs 
push the lower bounds on the heavy Higgs masses, we simplify the analysis, and use the bare masses 
instead. 

From the analysis of the tree-level mass matrices derived in this Section, Eqs. (31,34,35,38), and 
again working in the Higgs basis, a few interesting patterns can be extracted for the behaviour of 
the several mass eigenstates. Beginning with the lighter physical states (i = 2 for both charged and 
pseudoscalar states), the bare masses are very similar to what one has in the MSSM: 


m? < | cos 2/3 | Mz, 

(m 2) 2 « K) 2 + m 2 , 


(mf) 2 « K) 2 + My ,. (42) 


In the above, the approximations encode the fact that we have neglected the mixing between the several 
states (in each scalar, pseudoscalar and charged sector). Taking the mixing into account would lead to 
cumbersome expressions, involving all the states, which could only be numerically evaluated. For the 
remaining (heavier) states the physical masses are dominated to a large extent by the diagonal entries 
in the corresponding mass matrices, i.e., to 2 3 _ 66 . Although the general case is that the presence of 
the off-diagonal terms (m 2 ^ and bij) does not substantially modify the mass spectrum, the situation 
can be distinct if either one of the latter quantities becomes close to the values of the diagonal entries, 
namely 

m 2 j « 0.9 m 2 ;, or b z j « 0.35 m %. (43) 

If we are in the presence of the above regimes, then more mixing between the states is produced, and 
one is ultimately led to the appearance of tachyons. The presence of tachyons is also triggered by 
increasingly larger values of tan /3. 

Even though we will not address experimental issues, we will adopt the following naive bounds 
(mimicking the MSSM) for the bare masses of the lightest states [38]: 

m? > 75 GeV, m\ > 91 GeV, > 80 GeV, (44) 


where the first bound translates in taking tan/3 > 3. 
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4 Yukawa interaction Lagrangian 

In agreement with the superpotential introduced in Section 2, the Lagrangian for the interaction of 
Higgs with up- and clown-quarks can be written as 

^Yukawa = - E E 4 4 1? u' L + H.C. , (45) 

*=1,3,5 *=2,4,6 

where Ql,r (Y- l ' d ) are vectors (matrices) in flavour space, and the quarks appearing above are inter¬ 
action (rather than mass) eigenstates, a prime being used to emphasise the latter. This Lagrangian 
gives rise to the quark mass matrices and to scalar and pseudoscalar Higgs-quark-quark interactions. 


The mass terms for the quarks read 

£mass =-d' R M d d' L - u' R M u u' L + H.C. , (46) 

M d = J2 ™*E> MU = E WiY 'i ■ ( 47 ) 

i=l,3,5 2,4,6 

The resulting up- and down-type quark mass matrices can be diagonalised as 

v% M* Vf = diag(mf), q = u,d, (48) 

while the mass and interaction (primed) eigenstates are related by the following transformations 

v l Ql = 9l V r q' R = q R q = u,d. (49) 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix is defined as 

bcKM = V 2* V R . (50) 


From the above equations it is manifest that in the quark sector, mass matrices and Yukawa couplings 
are in general misaligned, the only exception occurring for Yukawa couplings which are proportional. 
This misalignment translates into the impossibility of diagonalising both matrices simultaneously, and 
is the source of the existence of tree-level FCNCs, as we will briefly discuss. 

Let us turn our attention to the Higgs-quark-quark interaction Lagrangian. In the unrotated basis, 
the latter reads 

^Yukawa = y= E [ 4 Y i 4 + 4 Y d d' L + H.C.] 

^ ^ i=l,3,5 

~7= E [ 4 Y i U U ' L + 4 * 4 Y i U 'l + H.C.] , (51) 

v ^ i=2,4,6 

while moving to the quark and Higgs mass eigenstate basis the couplings are 

^Yukawa = - 4 E [ ( V <0 ab % 4 4 + * 0^ h? 4 d\ + H.C. 

v i=l,3,5 

- 4 E ( V «)o6 h j 4 4 + i ( w «)ab 4 4 4 + H.C. . (52) 

v ^ 2 = 2 ,4,6 
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In the above, a, b denote quark flavours, while i,j = 1,..., 6 are Higgs indices, with s (p) denoting 
scalar (pseudoscalar) mass eigenstates. The latter are related to the original states as h s = S a a, 
h p = S v ip, as from Eqs. (40,41). The scalar (pseudoscalar) couplings V (W) are defined as 

(V,)ii ={Slh (y'YfVpU, 

(V*Y?V?U, (53) 

with i = 1,3,5 (2,4,6) for q = d(u) and j = 1,...,6. The several rotation matrices appearing in 
the previous equation were defined in Eqs. (40,41,48). As a final remark, it is important to stress 
that the matrices Vl.r which diagonalise the quark mass matrices do not, in general, diagonalise the 
corresponding Yukawa couplings. Hence, both scalar and pseudoscalar Higgs-quark-quark interactions 
may exhibit a strong non-diagonality in flavour space. 

Finally, let us point out that in the Higgs basis R\ and A 2 do not have flavour violating interactions. 
This is straightforward if one recalls that the Higgs basis mimics an MSSM-extended model. In the 
basis of the Higgs physical states, all six (five) scalars (pseudoscalars) mix, and all play a role in 
mediating the FCNC processes. 


5 FCNCs at the tree-level 

In this section, we compute the tree-level observables (such as neutral meson mass differences and CP 
violation in neutral meson mixing) induced by the exchange of neutral Higgs. These effects are absent 
in the SM, and play a determinant role in constraining the free parameters of the model. We discuss 
these effects for the case of the neutral kaons, as well as for the Bd , B s and D° — D° systems. 


5.1 A-meson oscillations and contributions to A mx 


We begin with the computation of the contributions of the several Higgs fields to the mass difference 
of the long- and short-lived neutral kaon states. In terms of effective Hamiltonians, the neutral kaon 
mass difference is defined as 


A rriK = tuk l — ttiks — 2 A4^ 


2 


(K° H 


AS =2 

eff 



(54) 


where is the effective Hamiltonian governing AS = 2 transitions. The Hamiltonian can be 

decomposed as 

A e A /= 2 = A tree + Aloop • (55) 

In the above we have separated the contributions arising from tree-level diagrams from those associated 
with box and higher-loop diagrams. We will focus on the tree-level contributions to A mx induced by 
the exchange of scalar and pseudoscalar Higgs. From the interaction Lagrangians previously derived, 
it is now straightforward to compute the effective Hamiltonian for the diagrams in Fig. 6 . One thus 
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Figure 6: Feynman diagrams contributing to Aitik at tree-level. 


has 7-t e ff = H° s + with 


njCT 

H-eft — 
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16(toH 2 

j=l —6 J 


- E 

2 = 1 ,3,5 




v d y 2 i + (vj - v d yi 75 ' 


+ H.c., 


(56) 


■yjV . 
n eff 


j=2-6 


16(m^) 2 


i=l,3,5 




tW, 


H.c.. 


(57) 


where the scalar and pseudoscalar Higgs masses, m s ’ p , have been defined in Eq. (40). Therefore, the 
contribution to the kaon mass difference associated with the exchange of a scalar Higgs boson ( h s ) is 
given by 


E 1 


;=i-e v m P 2 


E (Vdii* + Vd: 


n 2 

2O (^|(^)N)I^ U 


*=1,3,5 


- 

2 ) 

E (h4*-k4 {) 

2 =1,3,5 

<j? > \(ry B d)(s'y 6 d)\K 0 )' 


while the exchange of a pseudoscalar state ( h p ) reads 


(58) 


M 


K 

12 




=- Y 1 

8 ,^ 6 K ) 2 


E (W 


2 



(iF° |(*d)(5d)| A'°) 




2 1 

+ 

E (W^+Wd %) 

(K° \(s l5 d)(s l5 d)\ K°) 


2 =1,3,5 

J 


(59) 


The theoretical prediction for the value of A rriK thus obtained should be compared with the experi¬ 
mental value of (Awif)exp = 3.49 x 10 -12 MeV. We stress here that in a SM/MSSM-like scenario, the 
Wd, Vd matrices entering in Eqs. (58, 59) would be diagonal in flavour space, and thus no tree-level 
FCNC would occur. This is clear from Eqs. (53), since the flavour content of W and V is proportional 
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nriK 

497.6 MeV 

[38] 

rri£) o 

1.864 GeV 

[38] 

A rriK 

3.49 x 10" 12 MeV 

[38] 

Am do 

< 46.07 x 10" 12 MeV 

[38] 

Ik 

159 MeV 

[39] 

Id 

224 MeV 

[40] 

m Bd 

5.279 GeV 

[38] 

WB, 

5.369 GeV 

[38] 

A ^B d 

3.304 x 10 -13 GeV 

[38] 

A m Bs 

> 94.8 x 10” 13 GeV 

[38] 

fs d 

215 MeV 

[41] 

fs s 

245 MeV 

[41] 


Table 1: Numerical values used throughout the computation with the corresponding references. 


to (Vft Yf V^) a b- In the SM and MSSM, the matrices that diagonalise the quark mass matrices Vl,r 
also diagonalise the Yukawa couplings, so (V£ Yf V^) a b oc S a b- In multi-Higgs doublet models, where 
the underlying theory implies that the Yukawas are proportional among themselves (i.e. Yf 1 oc V)’ 
we also encounter a situation where no tree-level FCNCs emerge. 

Given the purpose of the analysis, we have adopted a simple approach regarding the computa¬ 
tion of the meson matrix elements, using the vacuum insertion approximation with non-renormalised 
operators. Following Ref. [39], we have 


(K \(sd)(sd)\K°) 

(K° I(s7 5 d)(s 7 5d)| K°) 


1 

12 


_L ( m ' K 

12 \m s + m d 


2 ' 


m K , 


1 

12 


11 / rn K \ 2 

12 \m s + m d J 


m K fx ) 


(60) 


where ttik is the kaon mass, and fx the kaon decay constant. In Table 1, we present several relevant 
input parameters for the computation of meson observables. 

Before concluding the analysis of the neutral kaon sector, let us mention that in a scenario where 
one has flavour violating neutral Higgs couplings, and given the most generic possibility of having a 
complex CKM matrix, it is natural to expect the occurrence of indirect CP violation at the tree-level. 
Although we will not pursue this issue in the following numerical analysis, let us just stress that the 
contributions to £k (which parametrises indirect CP violation in the kaon sector) are given by 


e”/ 4 Im [M^ 2 A 2 
y/2 |A„| 2 A m K 


(61) 


where A„ is defined from CKM elements as X u = V* s V ud , with as computed in Eqs. (58,59) (under 
the assumption of complex Yukawa couplings). This new tree-level contribution would have to be 
compatible with the SM loop contribution and the experimental bound £k = (2.284±0.014) x 10 -3 [38]. 


5.2 Other neutral meson systems: B d , B s and D° 

Computing the mass differences of neutral B and D mesons introduces no new elements into the 
analysis. The approach is entirely identical to that of the kaon system, the only difference lying 
in replacing the K° ( sd ) constituent quarks of Fig. 6 by (tic), (bd) and (bs), for D°, B d and B s , 
respectively. In each case the effective Hamiltonian for scalar and pseudoscalar Higgs exchange reads 
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as in Eqs. (56, 57), providing that the following replacements are done 6 : 

B d : V and W indices (21) —> (31). 

B s : V and W indices (21) —► (32). 

D: V and W computed for the up-sector (see Eq. (53)); V and W indices (21) —» (12); sum over 
interaction eigenstates i = 1, 3, 5 —> * = 2,4,6. 

In each case, the hadronic matrix elements should be also recomputed, and the predictions for 
each of the above processes should be confronted with the experimental data summarised in Table 1. 
Before concluding this section, let us briefly comment on the several observables mentioned here. 
First of all, it is widely recognised that, in models with tree-level FCNC, the most stringent bounds 
are usually associated with Am/<-. Regarding the B d sector, it has also been argued that in the 
absence of a predictive theory for the Yukawa couplings, the bounds associated with Ams d should 
be considered as a more reliable constraint, since they do not involve the mixing between the first 
two generations [19]. Although we include it in our analysis, the B s mass difference is not expected 
to add any new information. In the SM, this mixing is already maximal, and the addition of a new 
contribution would have little effect, the only exception occurring if new contributions matched exactly 
those of the SM, but had opposite sign, in which case a cancellation could take place. Still, this is a very 
fine-tuned scenario, and hardly significant, given the uncertainties associated with the computation. 
Finally, we turn our attention to the D° mass difference. As pointed out in [10,13] and [42], models 
allowing for FCNC at the tree-level may present the possibility of very large contributions to Amo, 
thus emerging as an excellent probe of new physics effects. The latter contributions are often harder 
to control than those associated with, for example, A uik- In fact, and as discussed in [13], the 
contribution of tree-level FCNCs to Am o could even exceed by a factor 20 those to A ttik- On the 
other hand, mixing in the D° sector is very sensitive to the hadronic model used to estimate the 
transition amplitudes, and there is still a very large uncertainty in deriving its decay constants, etc. 
Therefore, the constraints on a given model arising from Amo should not be over-emphasised. 

Another interesting issue is that of rare decays. It has been argued that, again when no theory for 
the full Yukawas is available, some rare decays become very sensitive to flavour changing contributions 
induced by Higgs exchange at the tree-level. In Ref. [19], the authors have identified that the most 
promising decay modes involve the leptonic sector, together with the H-mesons, and are /i —> ey, 
B d — > A'/.tr and B a — ■> /it. We will not pursue this issue in the present work, reserving some of the 
above modes (namely /i —■» ey) for a forthcoming analysis [43]. 

6 Flavour violation in six Higgs doublet models: results and 
discussion 

6.1 Yukawa couplings: the simple Fritzsch scheme 

We have carried out in the previous Section a general analysis of FCNC contributions. This can be 
applied to any particular model with three Higgs families, provided that the Yukawa couplings are 

6 In all cases, we are only computing an estimate value, not taking into account neither theoretical uncertainties (as 
those associated with the computation of the matrix elements), nor experimental errors. 
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known. Lead by simplicity, and following the analysis of [13], we take as an illustrative example the so- 
called “simple Fritzsch scheme” [31,32]. Taking this ansatz for each of the Yukawa couplings presents 
two main advantages: it leads to mass matrices with a fairly hierarchical structure, and allows to fit 
experimental data on quark masses and mixings for a reasonable number of free parameters. 


The “simple Fritzsch scheme” essentially consists in having all the weak quark eigenstates adopt 
identically structured couplings, which display suppressed flavour-changing elements for the first two 
generations in each family. More explicitly one has 



/ 0 

Ai 

V o 


Ai 0 \ 

0 Bi , 
Bi C * ) 


where i = 1,..., 6, and the entries are given by 7 


(62) 


Ai = at m\ ml, Bi = bi y ml ml > C* = c,: ml , (63) 

with a-i, bi and c,; coefficients of order one, and the quark masses of the j th generation and q = d (u) 
for i = 1, 3, 5 (2,4,6). One further advantage of this choice, which will become more evident in the 
subsequent discussion, is that under the above ansatz, the quark masses and mixings are independent 
of the chosen VEV regime (although the Yukawas are not). This means that for a given successful 
choice of the parameters a^, bi and Cj, one can study a number of distinct VEV schemes, corresponding 
to different tan/3 scenarios. 


Let us now present a specific numerical example. Taking the values for the input quark masses, 
mf entering in the ansatz of Eq. (63) as m u = 3 MeV, m d = 6.5 MeV, m c = 1.25 GeV, m s = 0.1 
GeV, rn t = 178 GeV, m b =4.5 GeV, and using the following choice of coefficients 


a = {0.45,0.30,0.35,0.45,0.20,0.50}, 
b = {0.25,0.40,0.45,0.30,0.20,0.35}, 

c = {0.30,0.45,0.30,0.30,0.40,0.25} , (64) 


one obtains the following set of output quark masses 


m u = 0.0042 GeV, m c = 1.37 GeV, rn t = 179.4 GeV 
m d = 0.0073 GeV, rri s = 0.087 GeV, rn h = 4.5 GeV 

and the associated mixing matrix 

/ 0.9742 0.2255 0.0031 

|Vckm| = 0.2251 0.9733 0.0438 

v 0.0129 0.0420 0.9990 


(65) 


( 66 ) 


which is in good agreement with experimental data. Having the relevant data (Yukawa couplings and 
CKM matrix), we can now proceed to estimate the contributions to neutral meson mass differences. 
As stressed before, we are considering a general multi-Higgs model, so that we have no definite scheme 
for the /x-terms and soft SUSY breaking masses. Therefore, and as done in Section 2, we will assume 
simple textures for the Higgs sector. 

7 For simplicity, and since our main concern is to illustrate the contributions to meson mass difference, we will assume 
that the Yukawa couplings are real, so that we will not examine the contributions to CP violation observables. 
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6.2 Tree-level FCNC in neutral mesons: numerical results 


The first step in evaluating the contributions to tree-level FCNC is to parameterise the Higgs sector, 
and thus obtain the mixing matrices and mass eigenstates entering in Eqs. (58,59). As discussed 
in Section 2, the Higgs basis is far more intuitive to work in than the original superpotential basis. 
Working in Higgs-basis notation (see Eqs. (14-24)), let us assume that the free parameters mfj and 
bij obey the following simple patterns: 


/ 0 

0 

® \ 

( ® 

0 

0 \ 

0 

*3 

V xlTeV, = 


*4 

y x ITeV, 

v ® 

y 

*5 / 

\ 0 

y 

*6 ) 


yfbTj = Vbx ITeV, (67) 


where the (g> denotes an entry which is fixed by the minimisation conditions (cf. Eqs. (22,23)). The 
above parametrisation can be further simplified by taking *3 = X 4 and *g = * 5 . Therefore, we have 
the following free parameters associated with the Higgs sector: 


X 3 , *5 , Vb, tan /3, 


( 68 ) 


and the pattern of VEVs, namely whether they are degenerate or hierarchical. In the numerical 
examples we will thus consider the following textures for my and bij. 

(A) : *3 = 1 , £5 = 5 , y = 0.4, Vb = 0.3 , 

(B) : a: 3 = 10, *5 = 10, y = 0.6, Vb = 0.5 . (69) 

In each case, we will consider several values for tan/3, taking it to lie in the range 3 < tan/3 < 12. 

Moreover, for every tan/3 we take to very distinct schemes of VEVs, that while respecting the EW 

symmetry breaking conditions and being in agreement with the chosen tan /3 value, exhibit either a 
degenerate pattern (w\ = W3 = w$, 1*2 = W4 = Wq) or a strongly hierarchical one ( 1*5 = 10 1*3 = 
100 w\, wq = 101*4 = lOOu^)- 

For the above range of tan/3, the Higgs spectra (lightest and heaviest scalar, pseudoscalar and 
charged states) is within the following ranges 8 : 

Texture (A) : 

79 GeV < mf < 90 GeV, 198 GeV < m\ < 409 GeV, 417 GeV < mf < 203 GeV, 

5009 GeV < m\, mg, mf < 5011 GeV. 

Texture (B) : 

79GeV <m{ <90GeV, 908GeV < m\ < 1547GeV, 909 GeV < mf < 1549GeV, 

10043 GeV < m s 6 , mf, mf < 10047GeV. (70) 


Let us begin by addressing the observable which is associated to most stringent constraints for this 
class of models: A mx- In Fig. 7, we plot ratio Am^/(Ami<-)exp as a function of tan/3 for textures 
(A) and (B), considering in each case a pattern of degenerate (Deg) or non-degenerate (NDeg) VEVs. 
Fig. 7 clearly reflects the most problematic aspect of this class of multi-Higgs doublet models. Without 

®In both cases m | and the masses of the heaviest states grow monotonically with tan 6 . However, in the case of 
Texture (A), since the mixing induced by b and the diagonal mass X 3 are of comparable magnitude, as tan/3 grows, so 
does the mixing between the states, and thus, as pointed out in Section 3, the states h% and h.j can become lighter 
with increasing tan f). 
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Figure 7: Attik /(A?n_R-) e xp as a function of tan (3. Full (dashed) red line denotes texture (A) for 
degenerate (hierarchical) VEVs, while a black full (dashed) line refers to texture (B) with degener¬ 
ate (hierarchical) VEVs. For Texture (A) a cross marks the appearance of tachyonic states. The 
experimental bound is depicted via a dotted horizontal line. 

a symmetry forbidding some of the Yukawa couplings, and if the Yukawas themselves do not exhibit 
a strong hierarchical character, the contributions to the neutral kaon mass difference can only be 
brought down to the experimental value via a set of very heavy Higgses, as those of texture (B). A 
Higgs spectrum closer to the EW scale, with a typical mass scale of 500 GeV, would generate, for the 
case of degenerate (non-degenerate) VEVs, contributions to A nrix of around 30 (3 x 10 4 ) (Am^-) e xp- 
It is also manifest that smaller values of tan/3 favour smaller contributions. This is due to having 
the Yukawa couplings for the down quarks proportional to sec /?. The relevance of the VEV regime 
should also be emphasised, since the latter plays a very important role. Even though the Yukawas 
enter in the contribution to already rotated by the matrices that diagonalise the quark mass 

matrices, it is clear that the smaller the VEVs associated to the first and second quark generations, 
the more enhanced will be the (12) matrix elements. In the case of degenerate VEVs, all the Yukawas 
are identically suppressed/enhanced 9 . Naturally, assuming such a large scale for the soft-breaking 
terms potentially leads to a fine tuning problem. As pointed out in Section 2.3, and even though the 
discussion was conducted for a distinct basis, masses above the few TeV scale are in principle within 
range of a more than 1% fine tuning. 

When compared to some previous studies, our results are more severe. However, let us stress that 
in our analysis we have taken a few distinct steps. Firstly, and in comparison to the ansatze used 
in [13], our Yukawa couplings are quite different, since accommodating the current CKM matrix data 
leads to values of at , bi and Ci quite smaller than those previously considered. This implies that the 
present Yukawas are less hierarchic. Moreover, we have explicitly taken into account the values of 
the VEVs, and considered the contributions from the exchange of all physical scalar and pseudoscalar 

®This is also related to the specific ansatz for the Yukawa couplings. 
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Figure 8: AmBj/fAmsJexp as a function of tan (3. Line and colour codes as in Fig. 7. 

Higgses (not weak interaction eigenstates), taking into account not only their masses, but also their 
mixings. It is also important to mention that the values of the hadronic matrix elements have been 
revised in the past years. 

In Fig. 8, we present the contributions to the mass difference of the Bd mesons. As expected, in 
this case it is far easier to comply with the experimental bounds. For the case of degenerate VEVs, 
even the “lightest” texture (A) succeeds in complying with the experimental bounds throughout the 
whole range of tan/3 considered, while for the heavier Higgs set (B), with non-degenerate VEVs. 
compatibility is obtained for the low tan/I regime (tan/3 < 6). 

To complete our study, we display in Fig. 9 the same analysis for the case of the B s meson system. 
In this case, the experimental bound is a lower (rather than upper bound). As discussed in Section 5, 
the SM already has maximal mixing in the B s system, and even though the new contributions may 
be of identical magnitude, there is no reason to expect a cancellation between SM and tree-level 
contributions, which would be associated to an extremely fine tuned scenario for all the parameters 
involved. 

From the meson systems so far discussed, the most severe constraints arise, as expected, from the 
Kl — Kg mass difference. As a final remark, and following the discussion of Section 5, we briefly 
comment on the predicted scenario involving the D° mesons. A study similar to those conducted for 
the K° and B sectors is of less interest, since there is little dependence of Amn /(A?7i£)) eX p on tan/3. 
This is due to having the up-type quarks Yukawa couplings associated to csc/3, as seen from Eq. (62). 

Rather than presenting a plot, we will briefly summarise the situation. With the exception of 
Texture (A), with non-degenerate VEVs, which induces contributions to A mo around 10 times its 
experimental value, all other cases predict contributions below the experimental bound, ranging from 
10~ 3 to 10" 1 (Amojexp. 

To summarise, ansatze for the Yukawa couplings of the “simple Fritzscli” type, when generalised 
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Figure 9: Am/j s /(ArnjjJexp as a function of tan/3. Line and colour codes as in Fig. 7. 

to multi-Higgs doublet models, typically induce tree-level FCNC’s, and require a Higgs spectrum at 
least of order 10 TeV, in order to ensure compatibility with experiment. Notice however that the 
results presented here are very dependent on the assumed scheme for the Yukawa couplings, and are 
only to be taken as an illustrative example, in the absence of a full theory for the Higgs-quark-quark 
interactions. Predictive models for Yukawa interactions of extended Higgs sector are likely to account 
for distinct, and hopefully less severe, bounds [7]. 


7 Conclusions 

Models that predict family replication in the Higgs sector offer a very aesthetic and phenomenologically 
interesting scenarios. Although there is abundant motivation for extending the Higgs sector (from 
string constructions, for example), in most cases the viability of these models is challenged by the 
occurrence of potentially dangerous FCNCs at the tree-level. 

We have analysed the most general form of the SUSY potential with three Higgs families, studying 
its minimisation, and deriving the tree-level expressions for the neutral (scalar and pseudoscalar) and 
charged Higgses mass matrices. 

The main goal of our work was to derive a model-independent evaluation of the tree-level contri¬ 
butions to neutral meson mass differences. We have computed the most general expression for the 
tree-level neutral Higgs mediated contributions to the mass difference of the neutral mesons. We 
took into account the exchange of all Higgs states, included the effects of mixing in the Higgs sector, 
and made no approximation with respect to dominant/sub-dominant contributions. This analysis is 
completely general, and can be applied to any given model with three Higgs families, independently 
of its Yukawa structure. As an example, in Section 6 we have assumed a simple ansatz for the Yukawa 
couplings (in analogy to what had previously been done), and have considered the contributions of 
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two distinct Higgs spectra to the K°, Bd, B s and D° mass differences, finding that the strongest 
bound - which as expected arises from A vtik - requires a spectrum of order 10 TeV. 

We again remark that the results for the Higgs masses are strongly dependent on the specific ansatz 
for the Yukawa couplings. Other ansatze, that account for a stronger hierarchy in the quark sector 
and still accommodate experimental data on quark masses and mixings, may generate quite smaller 
contributions to A mx, and thus require a lighter Higgs spectrum. This possibility is the subject of a 
forthcoming work [7], where it is shown that a Higgs spectra of order 1 TeV can be accommodated. On 
the other hand, it is also possible that while generating smaller contributions to Aijik, the different 
ansatze induce larger A ms, or Amo- 

In addition to the contributions to neutral meson mass difference, within the quark sector there 
are other processes that also deserve further investigation. For example, let us mention, at the one 
loop level, the very suppressed SM and MSSM B s decays. Processes involving the lepton sector also 
offer an even wider field for testing the new contributions induced by the additional Higgses (neutral 
and charged) [43]. CP violation, given the potential tree-level contributions to £k, may also become 
a stringent bound. 
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